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ABSTRACT 

According to the standard cosmological scenario, superclusters are objects that 
have just passed the turn around point and are collapsing. The dynamics of very few 
superclusters have been analysed up to now. In this paper we study the supercluster 
SC0028-0005, at redshift 0.22, identify the most prominent groups and/or clusters 
that make up the supercluster, and investigate the dynamic state of this structure. For 
the membership identification, we have used photometric and spectroscopic data from 
SDSS-DRIO, finding 6 main structures in a flat spatial distribution. We have also used 
a deep multi-band observation with MegaCam/CFHT to estimate de mass distribution 
through the weak-lensing effect. For the dynamical analysis, we have determined the 
relative distances along the line of sight within the supercluster using the Fundamental 
Plane of early-type galaxies. Finally, we have computed the peculiar velocities of each 
of the main structures. The 3D distribution suggests that SC0028-005 is indeed a 
collapsing supercluster, supporting the formation scenario of these structures. Using 
the spherical collapse model, we estimate that the mass within r = 10 Mpc should 
lie between 4 and 16 x IO^^Mq. The farthest detected members of the supercluster 
suggest that within ^ 60 Mpc the density contrast is <5 ^ 3 with respect to the critical 
density at z = 0.22, implying a total mass of ^ 4.6-16 x IO^^Mq, most of which in 
the form of low-mass galaxy groups or smaller substructures. 

Key -words: cosmology: theory - galaxies: supercluster. 


1 INTRODUCTION 

In the hierarchical paradigm of structure formation, the 
smallest structures are the first to collapse and virialize. 
They later get collected into progressively larger structures, 
each of which goes through the same stages of collapse and 
virialization. At present, the largest collapsed structures are 
clusters of galaxies, whereas superclusters are expected to 
be in the stage o f gravitational collapse, at least in their in¬ 
ner t e ns of Mpc (i Reisene gger et al]l2000l : iBatiste fc Batuskil 
l2013l : iMerluzzi et al.l 2015h . In the currently favoured cos¬ 
mological scenario, dominated by a cosmological constant 
or another form of “dark energy”, superclusters are start¬ 
ing to recede from each other at an accelerated rate, which 
will not allow them to get collected into even larger struc¬ 


tures. Therefore, they are the la rgest structures that wil l 
ever collaps e and virialize (e.g. Jn^ ;amnie&J;oey_l200^ 


Bushnet_^ 200dj_ Diinner et al.l l2006l . 20071 : Proust et alJ 
20061 : Rines et al. l2013h . If this scenario is correct, present 


superclusters play a pivotal role in our understanding of the 
evolution of the universe. Thus, it is very important to deter¬ 
mine their quantitative properties such as masses, sizes, and 
densities reliably, as well as testing this scenario as best as 
possible. Since superclusters have not yet collapsed and viri- 
alized, they do not stand out against the background den¬ 
sity field as clearly as individual galaxies or even clusters of 
galaxies (in terms of mass, optical light or X-ray emission), 
and equilibrium techniques used to determine the masses of 
clusters (hydrostatic equilibrium of the gas and virial equi¬ 
librium of the galaxies) do not apply to them. 


© 0000 RAS 



























2 O’Mill et. al. 


The study of superclusters has a long history, beginning 
with wo rks such as the identif ication of our own Virgo super¬ 
cluster (Ide Vaucouleurslll953t i and the statistical description 
of “clustering of sec ond order” (which ar e known nowadays 
as the superclusters: l^evman et al.lll95^ . Other studies ad¬ 
dressed the identifica tion, distribution and characterization 
of these objects fe. g. J^veer^ji^jjj^ 197^ Ta,go et al.lll984 
lEinasto et alJ Il984l . Il994 ; Zucca et alJ 1993li ! establishing 
that they are separated by large voids and clusters inside 
them are usually organized in chain-like structures. 

Some studies considered the morphology of superclus- 
ters, based on the spatia l distribution of their galaxies. 
lEinasto et al.l ll2007l . [^OlJ) analyzed a sample drawn from 
SDSS DR7, finding two main morphological types: filaments 
and others with a more complex, multibranch fine struc- 
ture. The wide mor phological variety of superclusters let 
lEinasto et al.l (1201 ill to suggest tha t their evolution have 
been dissimilar. Costa-Duarte et al.l (1201 il l used a kernel- 
based density field method to identify the superclusters and 
Minkowski Functionals to quantify their shape. They found 
that filaments and pancakes represent distinct morpholog¬ 
ical classes of superclusters in the Universe. The filaments 
tend to be richer, more luminous and larger than pancakes. 
It is then plausible to think that pancakes evolve towards 
filaments. 

While there is plenty of morphological analysis, dy¬ 
namical studies of galaxy superclusters are more scarce. 
In the hierarchical model of structure formation we expect 
that galaxy clusters form through accretion of galaxies and 
groups or merger with other clusters in the environment 
of superclusters. In the case of a supercluster with several 
clusters, it is possible that massive clusters grow through 
the gravitational collapse of the central parts of superclus¬ 
ters. Therefore, superclusters are still far from equilibrium 
today. The collapse scenario predicts that, among the clus¬ 
ters in a supercluster, those with higher observed redshifts 
are falling in to the centre from the front side (and thus 
are closer to us), while those with lower observed redshifts 
are falling in from the back (and are therefore more dis¬ 
tant). Consequently, in the gravitationally collapsing re¬ 
gion, the Hubble relation is reversed. This ef fect is indeed 
present in the core of the Shapley Supercluster (IProust et al.l 


20061: iRagone et ^l200d : iDiinner et akllioo^ . at z ~ 0.04. 
Batiste fc Batuskil 1120131 1 made a dynamical analysis of the 


Corona Borealis supercluster at « ~ 0.07 using data from 
the SDSS. They find this supercluster has broken from the 
Hubble flow and, based on dynamical simulations, conclude 
that a significant fraction of mass should reside outside the 
cluste rs comprising the supercluster. Recently, iTullv et al.l 
1 I 2 OI 4 II have used peculiar velocities to identify and describe 
a large and massive structure, Laniakea, comprising most of 
the galaxies in the local universe. 

The picture of gravitationally collapsing central parts of 
superclusters is required by the present cosmological model, 
but it has been tested only at low redshifts. Indeed, pre¬ 
vious dynamical studies, as those described above and in 
Sec. O dealt with superclusters below z = 0.1. In this work, 
we analyze the supercluster SC0028-0005, at 2: = 0.22, in 
order to test the collapse scenario for structures at inter¬ 
mediate redshift. This paper is organized as follows. In Sec¬ 
tion [ 2 ] we describe what is known about this supercluster 
and what is the data that we will analyze in this paper. 


In Section |3] we present an identification of the superclus¬ 
ter members. Section |4] contains a weak-lensing analysis of 
the substructures found in the supercluster. In Section [S] 
we show that the dynamical behaviour of the superclus¬ 
ter components is consistent with a collapsing scenario. For 
this analysis we use the Fundamental Plane of early-type 
galaxies to obtain relative distances and peculiar velocities 
of these components. Finally, our results are summarized in 
Sec. 15] Throughout this paper we adopt, when necessary, a 
ACDM cosmological model with Dm = 0.30, Da = 0.70 and 
Ho = 70 I 170 km s“^ Mpc“^. 


2 THE SUPERCLUSTER SC0028-0005 

In this study we focus on the galaxy supercluster SC0028- 
0005 (hereafter SC0028 for simplicity) at a = 00h28m and 
S = —00°05' (J2000.0). It is in the catalogue of iBasilakosI 
1 I 2 OO 3 II . obtained from the Sloan Digit al Sky Survey (SP SS) 
“Cut & Enhance” cluster catalogue (iGoto et al.l l2002fl by 
applying a percolation radius Rpc — 2Qh~^ Mpc. Among 
the 57 superclusters of this catalogue, SC0028 is the number 
38, with an assigned redshift z = 0.197, and is cl assified as 
a filam ent by using the shape finder estimator of IBasilakosI 
(I 2 OO 3 II . This structure was ch osen for thi s study because it 
was not too complex in the iGoto et al.l (I 2 OO 2 I 1 catalogue, 
containing only 3 galaxy clusters with relative distances of 
4 to 10h~^ Mpc between them, what was suggestive of the 
presence of strong peculiar motions. 


2.1 Spectroscopic and photometric samples 

We selected from the Data Release 10 of Sloan Digital Sky 
Survej0 (SDSS-DRIO) a sample of galaxies with spectro¬ 
scopic and photometric data within an area of radius 1.2 
degrees, centred on the supercluster. 

SDSS DRIO covers an additional 3100 sq. degree of sky 
over the previous release and includes spectra obtained with 
the new spectrographs APOGEE and BOSS, whose sky cov¬ 
erage includes the region of SG0028. 

The spectroscopic and photometric data sets were ex¬ 
tracted from the SpecObjAll and PhotoObj tables of the 
CasJob^ database. From the PhotoObj table we downloaded 
the objID, the coordinates, the model g and r magnitudes, 
the Galaxy extinction values, the effective radii derived from 
the de Vaucouleurs profile fit, as well as model magnitudes 
and axial ratios of the de Vaucouleurs fits, and the photo¬ 
metric redshifts. From the SpecObjAll table we downloaded 
the spectroscopic redshifts and the galaxy central velocity 
dispersions plus the respective errors. 

From this data we selected two galaxy samples. The first 
sample comprises 4921 galaxies with spectroscopic redshift. 
This sample is employed to identify the supercluster mem¬ 
bers and substructures, as well as in the dynamical analysis 
of Section [S] 

The second sample contains 35757 galaxies with photo¬ 
metric redshifts (including the galaxies of the previous spec¬ 
tral sample), selected from the PhotoObj table. This sample 


^ https://www.sdss3.org/drl0/ 

^ http://skyserver.sdss3.org/CasJobs/ 
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Figure 1. Redshift distribution within an area of 1.2 degrees radius centred on SC0028. Left PANEL: Spectroscopic redshift distribution. 
The inset shows the redshift distribution of galaxies in the Az < 0.02 interval around the supercluster mean redshift, as well as a Gaussian 
fit of this distribution. Right panel: Photometric redshifts. The gray histogram correspond to the redshift sample, reported here for 
easy comparison. In both panels, the dashed vertical lines correspond to the supercluster mean redshift and the magenta dotted lines 
show the galaxies within a redshift interval around the mean of Az < 0.02 for the SI sample and Azp/^^t < 0.04 tor the S2 sample. 


Table 1. CFHT Imaging characteristics 


Band 

Exposure 

(h) 

Seeing 

(arcsec) 

Completness 
(AB mag.) 

g 

1.5 

0.52 

24.5 

r 

3.3 

0.45 

24.6 

i 

2.5 

0.45 

24.3 


allowed an analysis of the photometric properties of the sub¬ 
structures identified spectroscopically and is also useful for 
their detection. 


2.2 CFHT Imaging 

The supercluster field was observed with CFHT/MegaCam 
in imaging mode, mostly for the weak-lensing analysis de¬ 
scribed in Sectionjd] We got deep, good quality g, r, i images, 
whose main features are shown in Table [T] 

The whole data processing (bias and overscan 
subtraction, flatfielding and sky subtraction), image 
combination and catalo gue extraction with SExtractor 
l|Bertin fc Arnoiit^ Il99fil i has been do ne by the Terapiijl 
team, using the procedures described in iBertin et al.l (l2002li . 

The total integrations were obtained with a series of 
~ 10 min. exposures per filter, with small position offsets, 
for the coverage of the cap between the camera CCDs. As 
a result, the final combined image in each band covers an 
area slightly smaller than one square degree (MegaCam’s 
full field of view). 


® http://terapix.iap.fr/ 


3 IDENTIFYING THE SUPERCLUSTER 

MEMBERS 

3.1 Selecting galaxies 

With the aim of selecting galaxies which are members of the 
supercluster, we start by defining Az as the absolute dif¬ 
ference between a given redshift and the supercluster mean 
redshift. For galaxies with spectroscopic redshifts we con¬ 
sider as members those with Az < 0.02. We call this sample 
SI. For those with photometric redshifts only, due to the 
larger uncertainties, we adopted Azphot < 0.04 (about twice 
the SDSS photometric redshift standard deviation). This we 
call the S2 sample. 

Left and right panels of Figure [T] show the resulting 
redshift distribution of the galaxies in the spectroscopic and 
photometric samples (Section 12.111 . respectively. The peak 
around z ~ 0.22 in these figures corresponds to the super¬ 
cluster. 

We selected galaxies from the spectroscopic sample 
which are supercluster members using the following ap¬ 
proach. We preselected all galaxies in the range 0.195 < z < 
0.245, a redshift interval large enough to not exclude possible 
members, and small enough to avoid most fore/background 
galaxies. We fitted a Gaussian to the resulting redshift dis¬ 
tribution, obtaining an average redshift of z = 0.220±0.001. 
Around this value we selected the spectroscopic sample, SI, 
using Az < 0.02, resulting in a total of 271 galaxies. For 
the photometric sample we adopted Azphot < 0.04, obtain¬ 
ing 3549 galaxies. In both panels of Figure [T] the magenta 
dotted lines bracket the membership ranges. 


3.2 Identifying substructures 

In order to identify substructures in our spectroscopic sam¬ 
ple, we present a density map of the SI sample in the upper 
left panel of Figure O This map was made by counting the 
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Figure 2. Substructures in the supercluster. Top-left: Surface density map of supercluster members in the SI (spectroscopic) sample. 
Top-right: Surface density map of supercluster members in the S2 (photometric) sample. Botton-left: Projected light distribution 
of the galaxies within photometric redshift Az < 0.04 of the supercluster mean spectroscopic redshift, z = 0.22. Bottom-right: Weak- 
lensing recovered mass distribution (See Section [4] for details). The green square corresponds to the 1 square degree MegaCam/CFHT 
field of view. In all panels, the red circles {R = 11.7 arcmin) correspond to maximum projected density regions (based on the SI sample), 
recentered at the location of the corresponding brightest galaxy in each maximum density box. The circles with solid lines are the 
substructures that will be considered in this study (see text for details). 


number of galaxies in equal size square cells with 15 arc- 
sec on a side. We found 9 well defined peaks in this density 
distribution which are most likely galaxy groups or clusters 
candidate members of SC0028. The brightest galaxy in each 
peak was taken as the group/cluster centre, and all galax¬ 
ies in the sample within a 2.5 Mpc radius (~ 11.7 ar¬ 
cmin) of this centre were considered as group/cluster mem¬ 
bers. In this panel, the circles show the region around each 
peak ascribed to a certain group or cluster. Table [2] presents 
the identification, coordinates, and mean redshift of each 
of these substructures. Notice that they may be groups or 
galaxy clusters, but hereafter we will call them substruc¬ 
tures. 

The upper right panel of Figure [2] is similar to the left 
panel, but uses the larger S2 sample (the photometric one) 
instead of the SI sample. 

The light map (Figure (2] lower left panel) was made by 
strongly smoothing the S2 sample of galaxies by a bidimen- 


sional Gaussian with intensity proportional to the galaxy 
luminosity and width proportional to their half-luminosity 
radius, a = 100i?5o. In this way, we can take into account 
the different apparent sizes of the galaxies. This luminosity 
field was then binned in a 2D grid with 15 x 15 arcsec^ cells, 
producing an image. Since the mean i?5o is ~ 1.1 arcsec, the 
typical galaxy luminosity was spread in a radius of ~ 7 pixel 
(Icr), which we found adequate to describe the substructures 
in the light distribution. This map is in qualitative agree¬ 
ment with the SI density map, as there are peaks in most 
regions previously identified (solid line circles in figure [2]). 

The bottom right panel of figure 0 shows the weak- 
lensing recovered projected mass distribution (see Sec¬ 
tion 0. The same structures identified in the previous maps 
are also present in this mass map, within the CFHT field of 
view (green square in the figure). This reassures that these 
substructures are not projection effects, but rather actual 
mass concentrations. Substructure 3 is the only exception. 
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Table 2. Structures within the supercluster: Identification, 
coordinates, mean redshift and red fraction. 


Id 

a (J2000.0) 

[h m s] 

S (J2000.0) 

[O . 

Z 

fraction of red 
galaxies [%] 

1 

00 28 44.28 

-000547.16 

0.218± 0.012 

78 

2 

00 28 59.35 

-1-002000.64 

0.222±0.015 

66 

3 

00 27 09.43 

-1-001413.95 

0.217±0.017 

55 

4 

00 26 48.78 

-001609.73 

0.226±0.017 

69 

5* 

00 30 21.65 

-002705.62 

0.225± 0.015 

61 

6* 

00 30 56.81 

-001007.15 

0.225± 0.017 

57 

7 

00 25 36.75 

-1-004701.09 

0.220±0.017 

38 

8 

00 28 03.93 

-010042.28 

0.217±0.015 

30 

9 

00 27 14.99 

-1-003246.67 

0.217± 0.015 

46 

SC 

00 28 00.00 

-00 05 00.00 

0.22040.001 

81** 


XJ 

o 



15.8 16.6 17.4 18.2 19.0 19.8 20.6 21.4 22.2 23.0 



Note*: There is some overlap between regions 5 and 6 (Figure 
[3. After checking the (photometric) redshift of the 6 galaxies in 
the overlap region we decided that this area should belong to the 
structure in region 5. 

Note**: The fraction of red galaxies was determined with the 
S2 sample. 


mag r 

Figure S. g — r colour versus r-band magnitude diagram for the 
galaxies in the photometric sample S2. The continuous line at 
(g — r) = 1.1 divides the sample in red and blue galaxies (see text 
for details). 


since it is detected in redshift space, upper left panel, but 
not clearly visible in any of the others. This probably means 
that substructure 3 is the least massive region of them all. 

In figure [3 we can see that areas 5 and 6 overlap. After 
analyzing the photometric redshifts of the 6 galaxies in the 
area in common, we decided to consider them as a part of 
group/cluster 5. 

Some of those substructures might be due to projec¬ 
tion effects and, to obtain a sample of more reliable physical 
structures, we have investigated their red galaxy content, 
since it is well know that the presence of a red sequence 
is conspicuous in large physical groups and cluster s (e.g. 
lQemleri[l974l : iDresslenllQSol : iPostrnan fc Gelleilll984l ) . 

We present in Figure [3] the {g — r) versus r colour- 
magnitude diagram for galaxies in sample S2. This figure 
shows a clear concentration of red galaxies. In order to make 
an inventory of red objects, we have divided our sample in 
two using (g — r) = 1.1 as a cutoff based on the typical col¬ 
ors of galaxies at redshift 2 = 0.2 as given bv iFukugita et aP 
(Il995l) : ellipticals (g — r) = 1.31, lenticulars (g — r) — 1.13, 
Sab galaxies {g — r) = 1.01. 

The fraction of red galaxies, by this definition, is shown 
in Table [2] for each identified structure. 

Based on this information we decided to remove sub¬ 
structures 7, 8 and 9 from our subsequent analysis, as they 
all have a fraction of red galaxies smaller than 50%. We sus¬ 
pect that these substructures may not be bona fide bound 
groups and/or clusters of galaxies. We thus kept, for the 
present study, 6 groups/clusters which we considered dense 
and red enough to be safely considered as substructures of 
SC0028. Note that structure 3 appears strong in the spec¬ 
troscopic redshift map (Figure [2]), while in the photometric 
redshift and light maps it appears diluted. We have consid¬ 
ered this structure significant, since more than 50% of its 
galaxies are red. 

Figure |4] presents results for the spectroscopic sample, 
SI. In the left panel we show the redshift distribution for 


each substructure. The dotted lines marks its mean redshift 
and the continuous line is the mean redshift of the super¬ 
cluster. 

On the middle panel we show the {g — r) x r colour- 
magnitude diagram of supercluster members in this sample. 
Dots of a given colour identify members of a specific struc¬ 
ture. On the rightmost panel we show the distribution of the 
(g — r) colour. 


4 WEAK-LENSING ANALYSIS 

In this section we present a weak-lensing analysis of SC0028, 
including sample selection and the PSF mapping and cor¬ 
rection. We describe the process of 2D mass reconstruction 
and fits to the shear profile of individual substructures. No¬ 
tice that the area used for the weak-lensing analysis (based 
on the CFHT observation) is significantly smaller than the 
area used for the spectroscopic and photometric analysis of 
the other sections, based on SDSS data (see Figure [2))- 

4.1 Background sample selection 

The selection of the background galaxies (i.e. those affected 
by the gravitational lens effect) has been done by identify¬ 
ing the regions on a (g — r) x (r — i) colour-colour diagram 
where they are abundant and that have the smallest con¬ 
tamination by foreground and supercluster galaxies. In this 
section we consider as members of the supercluster the 326 
galaxies with spectroscopic redshifts with redshifts between 
0.18 and 0.28. As a complement, we have also used pho¬ 
tometric redshifts in the same interval, from the CFHTLS 
deep field catalog (hereafter D; Coupon et al. 2009). 

In Figure [5] (left panel) we show a. {g — r) x {r — i) 
colour-colour diagram of the D galaxies (contours) 

with different colours that identify them as foreground 
(z < 0.18; blue), supercluster (0.18 ^ z ^ 0.28; green). 
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Figure 4. Results for sample SI. Left panel: spectroscopic red- 
shift distribution of each group candidate; dotted lines represent 
the mean redshift of the group and the continuous line is the 
mean redshift of the supercluster. Middle panel: The {g — r) xr 
colour-magnitude diagram. Each dot colour represents a struc¬ 
ture. Black dots are galaxies in the area not associated to any 
structure. Right panel: The {g — r) distribution. The dashed 
horizontal line at (g — r) = 1.1 on the last two panels indicates 
our adopted separation between red and blue galaxies. 

(A color version of this figure is available in the online journal) 


and background( 2 : > 0.28; red). With the aid of this in¬ 
formation we traced a triangle on this diagram, defined by 
vertices (-0.11,0.3), (0.58,1.3) and (0.45,-0.16), that selects 
the regions more heavily contaminated by non-background 
galaxies. All the galaxies outside this selected region are our 
primary background galaxy sample. Using the D catalogue 
we have estimated that the contamination of the background 
sample by field galaxies is only 0.6% within the magnitude 
range 17 < r < 24.6. In Figure [5] (right panel) we show the 
colour-colour diagram of the SC0028 field data, with the 
same triangular region overplotted, where the points repre¬ 
sent the spectroscopic redshifts. Within the same magnitude 
limits we have, on the CFHT image, a space density of ~ 20 
arcmin”^ candidate background galaxies. 



■V 

o 

d 


- 0.04 - 0.02 0.00 0.02 0.04 

El 

Figure 6. Values of components of ellipticity, ei and 62 , before 
(pink) and after (purple) the correction described in the text, 
where (ei) = 10.10 X 10“®, CTei = 7.32 X 10“®, (e2> = 2.14 X 10“® 
and (Tej = 5-21 x 10“®. 


maining parameters are the amplitudes (A) and a parameter 
related to the area of the base ellipse (ab). 

We selected 1944 stars based on the FWHM (~ 0.45 arc- 
sec) of their PSF in the range (18 < r < 22), where they 
have good signal-to-noise, are easily separable from galax¬ 
ies, and show no signs of saturation. From this sample we 
can estimate the PSF at any point of the image. We did 
that by smoothing the spatial variation of the relevant PSF 
parameters (ei, €2 and ab) with a Gaussian filter 


PSFpar = 


ZjWjPSFp, 


with Wi = exp 


\2a^ 


( 1 ) 


where PSFpar represents one of the relevant parameters at 
a given point of the image and di is the projected distance 
from the point to the star. The free parameter here is a. 
We found a best value of 110 arcsec by minimizing the vari¬ 
ance between predicted (interpolated) values with the actual 
values measured from the stars (See Figure |5J. 

Next, we applied Im2shape to the galaxies, tuning it 
to perform a PSF deconvolution using the predicted values 
for each galaxy position. We then excluded all the galaxies 
with a composed ellipticity error -|- > 0.45, ending 

up with a sample of 54187 galaxies (~ 15 gal. arcmin~^). 


4.2 Shape Measurements 

The observed shape of background galaxies is a combi¬ 
nation of their intrinsic ellipticity, the shear caused by 
the lens plus distortions due to the atmosphere and the 
telescope-binstrument optics (the PSF). The latter can be 
mapped through the observed shapes of stars (point sources 
in practice), and we can use it to deconvolve galaxy images. 

We measured galaxy and PSF shapes using the 
Im2shape software jBridle et al.ll2002 ). which estimates the 
ellipticity of astronomical objects by modeling them as a 
sum of Gaussians with elliptical bases. Stars are modeled as 
a single Gaussian whereas galaxies require a sum of two of 
those components, however keeping the centroid (xc, yc) and 
ellipticity (ei, £ 2 ) of both components as the same. The re- 


4.3 Mass Map 

We reconstruct the 2D density map fro m the shear infor¬ 
mation using the Le nsent2 algorithm (iBridle et al.lll9^ : 
[Marshall et al.ir2002l) . which creates a convergence (k) map 
through a maximum entropy method from the ellipticity 
components of each background galaxy ei and £2 and their 
respective uncertainties. To prevent overfitting, the resul¬ 
tant mass distribution is convolved with a Gaussian kernel. 
For that we choose an scale of 150 arcsec that yields a sig¬ 
nificant result, as estimated using the Bayesian evidence, 
without compromising too much the spatial resolution. 

The convergence is translated in to a physical mass den¬ 
sity by multiplying it with the critical lensing surface mass 
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Figure 5. Left panel: Colour-colour diagram identifying the regions that encompasses the galaxies of D catalog (contours) and galaxies. 
Supercluster galaxies (0.18 ^ ^ 0.28) are in green, foreground ones {z < 0.18) in blue, and background galaxies (z > 0.28) in red. The 

yellow triangle encloses the region where most foreground and supercluster galaxies lie. Right panel: Colour-colour diagram for galaxies 
in the field of the supercluster SC0028 with spectroscopic redshifts (points). 




density 

Ecwt = — -(2) 

47rGAA/ ^ ^ 

where Di, Da, and Djs are the angular diameter distances 
from the observer to the lens, from the observer to the 
source, and from the lens to the sonrce. 

Using the photometric redshift of the deep CFHTLS 
catal ogues we estimated a n average ratio (Da/Dis) of 1.39 
(See ICvpriano et al.l l2004l . for a description of the proce¬ 
dure). Therefore, we obtained Ecrit = 3.28 x 10^® Mq/Mpc^ 
= 0.685 g/crrd. The resultant mass distribution is shown in 
the lower right panel of Figure [21 

4.4 Individual Masses 


Table 3. Weak-lensing masses of SC0028 substructures with the 
CFHT imaging field. 


ID 

M 200 (lOi-^Mo) 

1 

0.65 ± 0.40 

2 

2.04 ±0.58 

4 

1.21 ±0.53 

5 

1.82 ±0.88 


galaxy and gt,i is the model prediction, aaii is an error re¬ 
lated to the intrinsic ellipticity of the galaxies, for which 
we measured a value of ~0.3, and am,i is associated to the 
measurement error of ellipticities. The latter is estimated as: 


We determined the masse s of each substructure in the su¬ 
percluster by fitting NFW (iNavarro et al.l[l996l . ll99^ 1 model 
predictions to their shear profiles up to a radius of 10 arcmin 
centered in each optically identified structure. Other values 
for radius have been tested, however the estimated masses 
remain unchanged within the statistical errors. The analyti¬ 
cal expressions for the radi al dependence of the shear of the 
NFW model can be seen in I Wright et al.l (I 2 OOOI I. 

The NFW profile can be completely defined by two pa¬ 
rameters: the mass inside a region with density 200 times 
the critical density, M 200 , and a concentration parameter c. 
The latter is poorly constrained by our data and thus we fix 
it at the value of c = 6.79, which is appropriate for IO^^Mq 
halos at 2 = 0.2 (jPrada et al.ll20lT . 

We obtain the masses by minimizing the standard sta¬ 
tistical misfit between data and model: 




(ct.i 9t,i) 


(3) 


where et,i is the measured tangential component of the el¬ 
lipticity with respect to the substructure centre for the i-th 


<Tm — 



(4) 


where and (Jej are the uncertainties in both components 
of the ellipticity estimated by Im2shape in the fitting pro¬ 
cess. 

In Table O we show the masses of the four identified 
substructures which are in the CFHT imaging field. The 
mass of region 3 is consistent with zero so it does not appear 
in this table. Region 6 does not appear either because its 
centre lies outside the CFHT image. 


5 DYNAMICAL ANALYSIS 

In this section we verify if the dynamical behaviour of the 
substructures identified in the supercluster are consistent 
with the collapse scenario described in the Introduction. For 
this, we first determine the relative distances and peculiar 
velocities of these substructures. We also examine our results 
with a simple dynamical model. 
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5.1 Scaling relations: the Fundamental Plane 


In this section we determine the Fundamental Plane of early- 
type galaxies to obtain, in the next section, estimates of the 
distances to the substructures identified in the supercluster. 

The Fundamental Plane (FP) is an empirical rela¬ 
tion between the central velocity dispersion (cro), the phys¬ 
ical effective radius (-Ro) and the mean surface bright¬ 
ness (ne) within the effective radius l|T erlevic h et al.l[l98ll : 
iDressler et al.l[l9^ : lDiorgovski &: Da^ 198^ 1. The FP can 
be represented as: 


logio(7?o) =alogio(cro)-f &logio(7e)-bc, (5) 

where le is the mean effective surface brightness in lin¬ 
ear unities, fie = —2.5log^Q(7e). 

If the FP is to be used as a distance estimator, it 
is convenient to ad opt a direct fit for the coefficients 
llBernardi et al.l[2003l l. because it minimizes the dispersion 
in the physical radius Ro- Another option is by applying or¬ 
thogonal fits, but this approach is more useful for the study 
of the global properties of elliptical galaxies or to constrain 
their underlying physics. Hence, we use here a direct fit for 
the calibration of the Fundamental Plane to use it as a dis¬ 
tance indicator^ _ _ _ 

Following iBernardi et al.l (l2003l l and iHvde fc Bernardil 
ll2009ll we first renormalized the effective radii from the 
SDSS data (rgdss) using the ratio of of the minor and major 
galaxy semi-axes (?ab) to account for the ellipticities of the 
galaxies in our sample. 

Thus, we adopted as our galaxy radius: 


r — rgdss^^AB ■ (6) 

This was d one to avoid a bias due to the distribution of 
ellipticities JBernardi et al.lboo,^ '). 

Since the SDSS uses a fixed fiber size, the fibers cover 
different galaxy physical areas at different distances. There¬ 
fore, we need to take this into consideration and correct the 
velocity dispersion for the spectroscopic galaxies in our sam¬ 
ple. 

The apert ure corrections f o r early type galaxies were 
calcul ated by I Jorgensen et al.l JlQQsIl and IWegner et al.l 
lll999lj as follows: 

___ /Tfiber , 0.04 ,r,\ 

— ^sdss( ) ? (7) 

r 

where r is the correcte d radius (see above) and rfiber is 1 
arcsec (lAhn et al.ll2012ll . 

Finally, the surface brightness in a circle of radius r is 
defined as: 


/ro = m-I-2.51ogio(27rr^) , (8) 


where m is the apparent magnitude corrected by extinc¬ 
tion and fc-correction. In this work. A:-correct ions were calcu- 
la ted using the public l y ava ilable software k-CORRECT v4.2 
of iBlanton fc R.oweisI 1 200it l. The effective radius Ro is in 
physical units of h^o kpc. 

In this subsection and in the next, we work with the SI 
sample. Since we are interested here in early-type galaxies, 
we consider only red galaxies, those with {g — r) > 1.1. 

We verified that the resulting sample is c onsistent with 
the m orphological classification of GalaxyZoo llLintott et al.l 
ll201lll i. available for 52 out of the 57 galaxies of our sam¬ 
ple. We used for this end the assigned probabilities of being 



a*loglO(aQ) + b*loglO(Ig) + c 


Figure 7. The Fundamental Plane relation. Black box symbols 
correspond to all red galaxies {{g — r) > 1.1) within Az < 0.02 
of the supercluster mean, whereas red symbols show galaxies be¬ 
longing to the supercluster substructures. The dotted and dashed 
lines are the best linear fittings to the black and red symbols, 
respectively. The fittings are almost identical, making the dis¬ 
tinction of the two lines difficult. The rms dispersions are 0.094 
and 0.079 for the black and red dots, respectively. 


ellipticals {P{E)). For the 52 galaxies we found that 49 have 
P{E) > 0.8, 2 P{E) > 0.71 and 1 has P{E) ~ 0.2. Addition¬ 
ally, the distribution of apparent axial-ratios is consistent 
with that of bona fide elliptical galaxies extracted from the 
Third Reference Catalogue of Bright Galaxies (RC3) by G. 
de Vaucouleurs, A. de Vac ouleurs, H.G. Corwin , R.J. Buta, 
P. Fouque, and G. Paturel ijCorwin et al.l l| 19941 )). We found 
~ 10% of our galaxies having (b/a) < 0.6, where a and b 
are the apparent major and minor axis, similar to the RC3 
sample. We conclude that our sample of galaxies should be 
less than 10% contaminated by non-elliptical galaxies, with 
no impact on the FP fitting. 

Figure [7] shows the best FP fit. The black box dots rep¬ 
resent the 123 red galaxies of the spectroscopic sample with 
reliable central velocity dispersion measurements and within 
Az < 0.02 of the mean supercluster redshift. The red dots 
represent 57 galaxies within the six substructures we identi¬ 
fied. The best fit parameters for the whole supercluster and 
for galaxies in the substructures are virtually indistinguish¬ 
able, with a = 1.035 ± 0.036, and b = —0.775 ± 0.031 and 
a rms scatter of 0.0941 and 0.0787 for the whole and sub¬ 
structure samples, respectively. These v alues agree well with 
those obtained bv ISaulder et al.l (l2013l) for the SDSS-DR8. 


5.2 Determination of distances and peculiar 
velocities 

In this Section, we first estimate relative distances to each 
substructure and, together with mean radial velocities, we 
estimate their peculiar radial velocities. Since we have on 
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Table 4. Results for the best fit of the the Fundamental Plane 
for the supercluster as a whole and the substructures. Column (1) 
shows the identification of the structure in consideration. Column 
(2) presents the Fundamental Plane coefficients (see text for more 
details), with errors calculated via a boostrapping procedure. Col¬ 
umn (3) is the rms dispersion in the fits for the supercluster as 
a whole and the substructures. Column (4) gives the number of 
galaxies used for the FP fitting. 


Sample id 


rms 

N 

sc 

a=1.035 ± 0.036 
b=-0.775 ± 0.031 
c=-9.105 ± 0.006 

0.0787 

57 

1 

c=-9.108 ± 0.015 

0.055 

11 

2 

c=-9.103 ± 0.012 

0.061 

10 

3 

c=-9.098 ± 0.014 

0.062 

8 

4 

c=-9.127 ± 0.014 

0.064 

9 

5 

c=-9.127 ± 0.012 

0.056 

10 

6 

c=-9.101 ± 0.017 

0.065 

9 


Table 5. Fundamental Plane mean redshifts and peculiar ve¬ 
locities for each substructure. Column (1) gives the identification 
of the substructures. Column (2) gives the cluster redshift de¬ 
termined from the Fundamental Plane and Column (3) gives the 
derived peculiar velocity. Column (4) gives the error in the dis¬ 
tances determined from the rms dispersion of the global fit and 
the number of objects in each substructure. 


Id 

Zfp 

Vp [km/s] 

err[%] 

1 

0.219 

-489.810 

4 

2 

0.217 

1025.577 

4 

3 

0.215 

269.444 

5 

4 

0.228 

-611.016 

5 

5 

0.228 

-815.231 

4 

6 

0.216 

2036.244 

4 


towards us, while a positive velocity indicates that it moves 
away from us. 

A summary of the main results of this section is given 
in Table O 


average only 10 galaxies in the SI sample per substructure 
(see Table 21 7 we do not use the FP to directly measure 
distances to the supercluster components; instead, the zero- 
point offset of each sustructure with respect to the overall fit 
can be used to find relative distances between them, keeping 
the coefficients a and b of Eq. fixed. The values of the 
fitting coefficients are shown in Table 0] for the supercluster 

as a whole an d for the substructure s. _ 

F ollowing IPearson et al.l (|2014l'l and iBatiste & Batus'S 
ll2013t) . we assume that the structures are, on average, at 
rest with respect to the CMB. Specifically, we do not assume 
any peculiar motion for the supercluster centroid. Individual 
substructures offsets are shown in Table ID 

The error of individual galaxy distances estimated from 
the FP is A = In(lO) x rms; since the dispersion of the FP 
relation is ~ 0.079 (Table |4ll, A is around 18%. in a gree¬ 
ment with the values obt ained bv IPearson et al] (l2014h and 
IBatiste fc Batuskil (l2013l i. For individual distances, the error 
percentage is reduced, assuming a Poissonian distribution, 
by A/y/N (where N is the number of galaxies in the sub¬ 
structure) and is tipically 6%. 

The comoving distances were calculated using the dif¬ 
ference of zer o points and converted to redshift by the ap¬ 
proximation (|Peeblesl|l993l l: 

D = 1^(1 - « 4283(1 - 0.2252) zhjf Mpc , (9) 

Hq 2 

which is appropriate at 2 ~ 0.22. We assume qo = —0.55. 
The peculiar velocities, Vp, were measured for each substruc¬ 
ture using the difference between the redshift obtained with 
the Fundamental Plane (zfp) and the mean spectroscopic 
redshift of the supercluster (2m): 


Vp = c 


f Zm - Zfp \ 
\ 1 -I- Zfp ) 


( 10 ) 


A negative velocity indicates that the substructure moves 


5.3 SC0028 as a collapsing structure 

Now we test whether or not the substructures we identi¬ 
fied are in a process of collapse towards the supercluster 
barycentre. The most identifiable feature of the process is 
a simple redshift space distortion. Substructures which are 
physically closer to the observer than the supercluster centre 
would feel its gravitational pull towards it, and thus will ap¬ 
pear as having positive peculiar line-of-sight velocities (i.e., 
Vp points away from the observer, adding to the receding 
velocity.). The opposite is expected in case of substructures 
which are physically farther than the supercluster barycen¬ 
tre. 

To test this hypothesis, we have used the FP distances 
and peculiar velocities determined in the previous subsec¬ 
tion. The average distance of SC0028 is ~ 1090.50 Mpc 
so, disregarding statistical uncertainties, substructures 2, 3 
and 6 would be in front of the supercluster and thus should 
have positive peculiar velocities, while substructures 1, 4 
and 5 would be behind and are supposed to have negative 
peculiar velocities. 

Assuming a Normal distribution of probabilities for the 
peculiar velocities and distances, it is easy to estimate the 
total probability that a particular cluster/group is in one 
of the two expected quadrants in the Dr x Vp space, i.e. 
the probability of having positive peculiar velocities and 
negative relative distances, plus the probability of having 
negative peculiar velocities and positive relative distances. 
Assuming independence between those quantities and their 
errors, we find the following probabilities for each substruc¬ 
ture be indeed in the expected quadrant where it is found: 
from substructure 1 to 6 we obtain, respectively, 59.4, 53.5, 
53.2, 65.2, 73.2, 65.3 %. Combining all these p-values using 
the Fisher method (iHedees and 01kinlll985h . we obtain a 
probability of 92.2%. This number can be interpreted as a 
moderate evidence for the detection of the infall of these 
substructures. 

Clearly most of the signal comes from substructures 4 
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Figure 8. Left panel: Distances relative to the supercluster center versus peculiar velocities for the six substructures, using the Funda¬ 
mental Plane method. Right panel: The 3-dimensional spacial distribution of SC0028 substructures also showing the LOS velocity vector. 
The length of the Vp vectors are in arbitrary units, proportional to the velocities derived using the Fundamental Plane. 


and 5, which, in particular, presents no ambiguity regarding 
being on the far side of the supercluster. The overall proba¬ 
bility is fairly insensitive to the exclusion of the less certain 
substructure 3. 

The left panel of Fig. [8] shows the phase space D-Vp 
for the 6 substructures of SC0028. The corresponding three 
dimensional structure is depicted on the right panel of Fig.|Sl 
where the present collapsing state of the supercluster can be 
appreciated. The errors in Vp were calculated through error 
propagation considering the error in the distance of each 
substructure. 


5.4 Spherical collapse model 


In order to estimate the supercluster mass, we employ the 
spherical collapse model in an expanding universe, that de- 
scri bes the free-fall evolution of a shell encl osing some mass 
M (iGunn fc Gott|[l972l : iLahav et al.lll99lh . Albeit simple, 
the spherical collapse model has been extensively studied 
and used to estimate the mass o f superclusters by several 
authors (e.g. DimrieretalJ_ 200^ Arava-Melo et aP l2009l : 


iBatiste fc Batus'5 20131: Pearson et al.l 20141. 


Here, we follow the time evolution of a spheri- 
cal shell by solving s imultaneously the motion equation 
llPeebles fc Ratr^l2C)03ll . 


dV GM 

_ + __H = 

and the background Friedman equation, 
dz 


dt 


+ (1 + z)H{z) — 0, 


( 11 ) 


( 12 ) 


with 


H{z) = Ho^J^a + ^m{1 + z), 

For each one of the six substructures, we have solved 



Figure 9. Free fall trajectories for the six substructures in 
SC0028 in phase-space. The inicial conditions are marked as large 
dots, with the corresponding numbers. 


the above differential equations using, for the inicial condi¬ 
tions, z(0) = 0.22, and r(0) and Ur(0) corresponding to each 
group as in Fig.jS] Notice that Wr(0) is always negative since 
the velocity vectors for all substructures point towards the 
centre of SC0028. 

After setting the inicial conditions and cosmological pa¬ 
rameters (the standard ACDM model), there still one un¬ 
known in the equations since we do not know a priori the 
mass of the supercluster. We have searched interactively for 
a solution with the present observed values for r(0) and Vr(0) 
and with the smallest density contrast, Ssc = ’p^^/pc, where 
Pc is the critical density at the supercluster redshift 2 = 0.22. 
This yields a lower limit estimate for the supercluster mass. 

We have integrated the motion equations forward and 
backward in time, so that we had the full trajectories from 
the turn-around point onward. Figure shows these trajec¬ 
tories in phase-space and Table [6] summarizes the numerical 
integration results. 
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Table 6. Initial conditions, r(0) and Dr-(O), and results of the 
integration of the spherical collapse equation. 6sc is the mean 
density contrast of the supercluster with respect to the critical 
density. fxA and rxA are the time of turn-around and the radial 
distance at that time. 


Obj. 

Id 

r(0) 

[Mpc] 

Vr{0) 

[km/s] 

<^SC 

M[< r(0)] 
[IOISMq] 

Ita 

[Gyr] 

rxA 

[Mpc] 

1 

10.0 

-490 

6.1 

4.3 

-5.7 

11.1 

2 

5.5 

-1026 

21.6 

2.6 

-6.5 

8.0 

3 

14.1 

-269 

3.6 

7.2 

-5.5 

14.6 

4 

61.0 

-611 

2.8 

445 

-3.6 

61.9 

5 

60.0 

-815 

3.2 

489 

-5.7 

61.0 

6 

9.8 

-2036 

24.9 

16.4 

-7.4 

14.7 


Figure [5] shows that substructures 1, 2, 3 and 6 have 
nearby trajectories in phase space, while substructures 4 and 
5 have very similar trajectories, but distinct from the other 
four. The latter pair seems to have about the same turn¬ 
around radius, compatible with an enclosing total mass of 
~ 4.5-4. 9 X Mq. This large mass, however, corresponds 
to a mean density contrast of only ~ 3 above the critical 
background density ai z = 0.22. 

The case of substructures 1, 2, 3 and 6 is more com¬ 
plex. Objects 2 and 6 suggest a high density contrast, ~ 22- 
25, while objects 1 and 3 suggest a low density contrast, in 
agreement with substructures 4 and 5 (between 5sc ~ 4 and 
6). On the other hand, substructures 3 and 6 have almost 
the same turn-around radius, rxA ~ 14 Mpc. 

Given the uncertainties in the distance determination 
and our lack of information on the transverse component 
of the substructures peculiar velocities - we have used the 
line-of-sight peculiar velocity as a rough estimate of the ra¬ 
dial peculiar velocity, Vr - we can only have an order of 
magnitude estimate on the supercluster mass. At a radius 
of 10 Mpc, the free-fall motion suggests that mass should lie 
between ~ 4 and 16 x IO^^Mq. 

The weak lensing masses of the individual substructures 
are of the order of 1-2 x (see Table. [3]). Within 

10 Mpc from the supercluster center, there are 3 to 4 of such 
substructures that adds up to about ~ 3 x 10^^ (without 
substructure 6, which lies outside the footprint of CFHT). 
Therefore, most of this supercluster mass budget should be 
in the form of lower than ~ 5 x lO^^M© galaxy groups. That 
is below our detectable level with the present data. 


6 SUMMARY AND CONCLUSIONS 

In this work we have analyzed the gravitational dynamics 
of the supercluster SC0028-005 with spectroscopic and pho¬ 
tometric data from SDSS-DRIO. This su percluster was se - 
lected from the supercluster catalogue of iBasilakosI (|2003h . 
where it was classified as a filamentary structure. 

We have chosen an area of 1.2 square degrees around 
the supercluster and proceeded to identify substructures us¬ 
ing density maps and red-sequence analysis. We identified 9 
substructures throughout the area, but 3 of them were dis¬ 
carded because their fraction of red galaxies was less than 
50%. 

To test the collapse scenario for supercluster evolution. 


we have determined the distances and peculiar velocities 
for each substructure in the supercluster. In order to study 
the dynamics of these substructures, the Fundamental Plane 
was used as a distance indicator. For the spectroscopic sam¬ 
ple, each substructure has, on average, 10 galaxies. This 
number is too small to determine the local Fundamental 
Plane parameters. Instead we have used the globally fitted 
slopes and then estimated zero point offsets for each of the 
six substructures. 

In order to obtain the individual distances and pecu¬ 
liar velocities, we considered that the supercluster’s centroid 
does not have peculiar movements and, knowing the super¬ 
cluster’s distance beforehand, we computed the individual 
distances for each structure. 

The resulting picture for SC0028 is of three substruc¬ 
tures approaching the supercluster barycentre from the far 
side and three from the near side, based on the computed 
empirical relations with photo-spectroscopic data. The spa¬ 
tial distribution and peculiar velocities of the detected mem¬ 
bers of SC0028 support the collapsing supercluster scenario. 
Assuming that the dynamics may be described by a simple 
free-fall spherical collapse model, the overall density con¬ 
trast is (5 ~ 3 inside a radius of 60 Mpc. In the inner region, 
inside 10 Mpc, the mass should be 4 — 16 x 10^® M0. Most 
of this should be in the form of low mass (M < 5 x 10^® Mq) 
substructures. 


7 ACKNOWLEDGMENTS 

We acknowledge financial support form CAPES/Cofecub 
French-Brazilian cooperation, CNPq and FAPESP. 

Funding for SDSS-III has been provided by the Al¬ 
fred P. Sloan Foundation, the Participating Institutions, 
the National Science Foundation, and the U.S. Depart¬ 
ment of Energy Office of Science. The SDSS-III web site 
is http://www.sdss3.org/ 

SDSS-III is managed by the Astrophysical Research 
Consortium for the Participating Institutions of the SDSS- 
III Collaboration including the University of Arizona, the 
Brazilian Participation Group, Brookhaven National Lab¬ 
oratory, Carnegie Mellon University, University of Florida, 
the French Participation Group, the German Participation 
Group, Harvard University, the Institute de Astrofisica de 
Canarias, the Michigan State/Notre Dame/JINA Partici¬ 
pation Group, Johns Hopkins University, Lawrence Berke¬ 
ley National Laboratory, Max Planck Institute for Astro¬ 
physics, Max Planck Institute for Extraterrestrial Physics, 
New Mexico State University, New York University, Ohio 
State University, Pennsylvania State University, University 
of Portsmouth, Princeton University, the Spanish Partic¬ 
ipation Group, University of Tokyo, University of Utah, 
Vanderbilt University, University of Virginia, University of 
Washington, and Yale University. 

This work is based (in part) on data products pro¬ 
duced at the TERAPIX data center located at the Institut 
d’Astrophysique de Paris and on observations made with 
ESO Telescopes at the La Silla Paranal Observatory (Chile) 
under programs ID 077.A-0262 and 085.A-0116. 


© 0000 RAS, MNRAS OOO.fTHT^ 











12 O’Mill et. al. 


REFERENCES 

Ahn, C. P., Alexandroff, R., Allende Prieto, C., et al., 2012, 
ApJS, 203, 21 

Ahn, C. P., Alexandroff, R., Allende Prieto, C., et al., 2014, 
ApJS, 211, 17 

Araya-Melo P.A., Reisenberg A., Meza A., van de Weygaert 
R., Diinner R., Quintana H., 2009, MNRAS399, 97 
Basilakos, S. 2003, MNRAS, 344, 602 
Batiste M., Batuski D.J., 2013, MNRAS, 436, 3331 
Bernardi, M., Sheth, R. K., Annis, J., et al. 2003, AJ, 125, 
1882 

Berlin, E., Arnouts S., 1996, A&AS, 117, 393 
Blanton, M. R., & Roweis, S. 2007, AJ, 133, 734 
Berlin E., Mellier Y., Radovich M., Missonnier G., Didelon 
R, Morin B., 2002, ASPC, 281, 228 
Bridle, S. L., Hobson, M. P., Lasenby, A. N. & Saunders, 
R. 1998, MNRAS, 299, 895 

Bridle, S., Kneib, J.-P., Bardeau, S.,Gull, S.F., 2002, in 
’The shapes of Galaxies and their Dark Halos’ Yale Cos¬ 
mology workshop, 28-30 May 2001, World Scientific. 
Busha, M. T., Adams, F. C., Wechsler, R. H., & Evrard, 
A. E. 2003, ApJ, 596, 713 

Costa-Duarte, M. V., Sodre, L., Jr., & Durret, F. 2011, 
MNRAS, 411, 1716 

Coupon, J., et al.,2009, A& A, 500, 981 
Corwin, H. G., Jr., Buta, R. J., & de Vaucouleurs, G. 1994, 
AJ, 108, 2128 

Cypriano, E. S.,Sodre, L. J., Kneib, J.-P., & Campusano, 
L. E. 2004, ApJ, 613, 95 
de Vaucouleurs G., 1953, AJ 58, 30 
Djorgovski, S., & Davis, M. 1987, ApJ, 313, 59 
Dressier, A., Lynden-Bell, D., Burstein, D., et al. 1987, 
ApJ, 313, 42 

Diinner, R., Araya, P. A., Meza, A., Reisenegger, A., 2006, 
MNRAS366, 803 

Diinner, R., Reisenegger, A., Meza, A., Araya, P. A., & 
Quintana, H. 2007, MNRAS, 376, 1577 
Dressier, A. 1980, ApJ, 236, 351 

Einasto, J., Klypin, A. A., Saar, E., & Shandarin, S. F. 
1984, MNRAS, 206, 529 

Einasto, M., Einasto, J., Tago, E., Dalton, G. B., & An- 
dernach, H. 1994, MNRAS, 269, 301 
Einasto, M., Saar, E., Liivamagi, L. J., et al. 2007, A& A, 
476, 697 

Einasto, M., Liivamagi, L. J., Tago, E., et al. 2011, A& A, 
532, A5 

Faber, S. M., & Jackson, R. E. 1976, ApJ, 204, 668 
Fukugita M., Shimasaku K., Ichikawa T., 1995, PASP 107, 
945 

Goto, T., Sekiguchi, M., Nichol, R. C., et al. 2002, AJ, 123, 
1807 

Gunn J.E., Gott J.R., HI, 1972, ApJ176, 1 
Hedges, L., Olkin, I. ,1985, Statistical Methods for Meta- 
Analysis, Academic Press, San Diego. 

Hyde, J. B., & Bernardi, M. 2009, MNRAS, 396, 1171 
Jprgensen, I., Franx, M., & Kjaergaard, P. 1995, MNRAS, 
276, 1341 

Joeveer, M., Einasto, J., & Tago, E. 1978, MNRAS, 185, 
357 

Kormendy, J. 1977, ApJ, 218, 333 

Lahav O., Lilje P.B., Primack J.R., Rees M.J., 1991, MN- 


RAS251, 128 

Lintott, C., Schawinski, K., Bamford, S., et al. 2011, MN¬ 
RAS, 410, 166 

Marshall, P. J., Hobson, M. R, Gull, S.F., Bridle, S. L. 
2002, MNRAS, 335, 1037 

Merluzzi, P., Busarelo, G., Haines, C. P., et al., 2014, MN¬ 
RAS, 446, 803 

Nagamine, K., & Loeb, A. 2003, New Astronomy, 8, 439 
Navarro, J. F., Frenk, C. S., & White, S. D. M. 1996, ApJ, 
462, 563 

Navarro, J. F., Frenk, C. S., & White, S. D. M. 1997, ApJ, 
490, 493 

Neyman J., Scott E. L., Shane C.D., 1956, in:“Proceedings 
of the Third Berkeley Symp. on Mathematical Statistics 
and Probability”, 3, 75 
Oemler, A., Jr. 1974, ApJ, 194, 1 

Pearson D.W., Batiste M., Batuski D.J., 2014, MNRAS441, 
1601 

Peebles P.J.E., 1993, “Principles of Physical Cosmology”, 
Princeton University Press. 

Peebles P.J.E., Ratra B., 2003, Rev. Modern Phys. 75, 559 
Postman, M., & Geller, M. J. 1984, ApJ, 281, 95 
Prada, F., Klypin, A.A., Cuesta, A.J., Primark, J. 2012, 
MNRAS, 423, 3018 

Proust, D., Quintana, H., Carrasco, E. R., et al. 2006, The 
Messenger, 124, 30 

Ragone, C. J., Muriel, H., Proust, D., Reisenegger, A., & 
Quintana, H. 2006, A& A, 445, 819 
Rines, K., Geller, M. J., Diaferio, A., Kurtz, M. J. 2013, 
ApJ, 767, 15 

Reisenegger, A., Quintana, H., Carrasco, E. R., & Maze, J. 
2000, AJ, 120, 523 

Saulder C., Mieske S., Zeilinger W., Chilingarian I., 2013, 
A& A, 557, id.A21 

Tago, E., Einasto, J., & Saar, E. 1984, MNRAS, 206, 559 
Terlevich, R., Davies, R. L., Faber, S. M., & Burstein, D. 
1981, MNRAS, 196, 381 

Tully, R. B., Courtois, H., Hoffman, Y., Pomarede, D. 
2014, Nature, 513, 71 

Wegner, G., Colless, M., Saglia, R. R, et al. 1999, MNRAS, 
305, 259 

Wright, C., O. & Brainerd, T., G. 2000, ApJ, 534, 34. 
Zucca, E., Zamorani, G., Scaramella, R., & Vettolani, G. 
1993, ApJ, 407, 470 


© 0000 RAS, MNRAS 000. [TIO 


